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We have studied the ground state of a simple cubic magnetic cluster, which contains a spin s 
at each corner site. The ground state of such cluster depends on the competition between nearest, 
next-nearest and next-next-nearest-neighbor exchange interactions. We have calculated "phase" 
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O , diagrams for the Izing clusters with s = — ; 1; — ; — ; oo and for the Heisenberg cluster with s = — ; 1. 
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QWe have found that the "phase" diagram is remarkably independent on the s value. It is important 
also that the Izing "phase" diagram can be used as a rough approximation for the Heisenberg model. 
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In recent years much attention has been given to mag- 
' netic properties of molecular magnets - clusters, each 
^ , containing relatively small number of paramagnetic ions. 
The best known examples of such clusters are Mni2 and 
Fes, which have been intensively studied|-due to quan- 
tum tunneling magnetization phenomenal!! The number 
of transition metal ions in molecular paramagnetic clus- 
i ters extends from 2 Jjn dimcrs) up to 30, as in polyox- 
ometalate Mo72Fe3oJ3 A spin quantum number s varies 
I , from 1/2, e.g. for polyolate-bridged copper clustersB, up 
to 5/2 for iron and manganese clusters. 

It appears that the majority molecular magnets can be 
described by a model of single-particle localized magnetic 
moments ("spins") coupled by the Heisenberg exchange 
interactions .□ In spite of smallness of molecular magnets, 
a numerical calculation of the energy spectrum could be a 
nontrivial task, since the dimension of the Hilbert space 

■ for a system of N spins s, given by (2s + 1)^, grows 
■<^J- \ rapidly with N and s. For example, in the case of Mni2 
CN| ■ the exact calculation of energies of all spin states arising 

from the coupling of eight s a = 2 spins and four Sb = 3/2, 
whose number is (2s a + l) 8 (2sft + 1) 8 = 10 s , is practicalbz 
out side the capability of the most powerful computers. u 
In magnetic clusters each of spin has a few neighbors, 
resulting in a frustration of exchange interactions. Be- 
cause of complexity of exchange pathways, both a sign 

■ and a magnitude of exchange constants are difficult to 
I predict a priori. A cluster becomes magnetic, if it has 

a non-zero-spin ground state. It is interesting to study 
effects of competing interactions on the ground state of a 
relatively small and highly symmetrical spin cluster, for 
which a wide range of exchange parameters can be rel- 
atively easily examined. In this paper we continue our 

■ examination of zero-temperature magnetic properties of 
cubic clusterso and present results of an investigation of 
ground-state properties of a cubic cluster of eight spins 
s, coupled by Izing or Heisenberg interactions. There 
are enough examples of real molecular clusters, contain- 
ing metal ions at corner sites of a simple cube. Un- 
fortunetly, these ions often have "non-magnetic" valence 
(e.g. copper-dithiolato species and related compoundsQ, 
cobalt cyanide clusters^). However, many potentially 
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interesting cubic clusters, (for instance, lanthanide sul- 
fido clustersu, met-carscj) remain unexplored by mag- 
netic techniques. 

A sketch of the cluster is shown in Fig. 1. The model 
Heisenbcrg Hamiltonian of the system takes the form 

H = J\H nn + JlH nmx + J3H nnnn (1) 

with 

H nn — SqSi + S0S3 + S0S4 + S1S2 + S1S5 + 5253 + 

-\-S2SQ + S3S7 + S4S5 + + S§Sq + SqSj; (lo) 
H nnn = S0S2 + S0S5 + S0S7 + S1S3 + S1S4 + 5i56+ 

+ S2S5 + S2S7 + ^3^4 + 53^6 + S^S^ + S5S7; (lb) 
Hnnnn = SoSq + S1S7 + S2S4 + S3S5; (lc) 

for the Heisenberg Hamiltonian, and 
H nn = SqS^ + 5g5f + 5q5| + 5f5f + 5f 5f + 5*| 5*1+ 

-\-S2SQ + 5*3 + 5^5^ + S^Sj + S^Sq + SqSy', (1^0 

QZ QZ I QZ QZ , QZ QZ , QZ QZ ■ QZ QZ ■ QZ QZ i 

-\-S2S5 + S 2 S^ + 5354 + S^Sq + 5*4 5g — (— *S* 5 iS* ^ ^ (lc) 

Hnnnn — 5p5| + S^S^ + S 2 S^ + 5|5g; (1/) 

for the Izing Hamiltonian. Here 5, is the spin operator 
located at lattice site i, and Sf denotes its z component. 
Terms J\H nn , J2H nnn and J-sHnnnn describes the NN, 
NNN and NNNN interactions, correspondingly. It is im- 
portant that the Heisenberg Hamiltonian (1) commutes 
with the operator of the total spin 5 2 as well as its z com- 
ponent 5 Z , where 5* = S0+S1 + S2+S3 + S4+S5 + S6+S7. 
Therefore, the energy level scheme consists of a number 
multiplets each of which can be described by quantum 
numbers 5 and S z ■ The number of states is easily derived 
by noting that the spin states of Si form a basis for the 
irreproducible representation D$- For example, in case 
of ,s = 1/2 the 256 states must form a basis for the repre- 
sentation D 1/2 xD 1/2 x D 1/2 x D 1/2 xD 1/2 xD 1/2 xD 1/2 x 
D1/2 = (Do + Dx) x (D0 + D1) X (D + D 1 ) x (D„ + I>i) - 
14L> + 28L>i + 20L» 2 + 7D 3 + L> 4 . Thus we expect 14 
singlet 5 = states, 28 triplets with 5 = 1 and so on. 
The eigenfunctions of (1) can be expressed as linear com- 
binations of basis functions ip m = |... ... [k ...}, 

7 

where m — ^ a^ 1 , = 1 or according as the i- 

i=0 

th spin is up (|i) or down (ji). All the eigenvalues e ra 
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and the corresponding eigenfunctions <f> n of the Hamilto- 
nian (1) can be found by the exact diagonalization of a 
256x256 matrix. (for example, using the standard House- 
holder methocO). It should be noted that the functions 
(p m are exactly equal to the eigenfunctions of the Izing 
Hamiltonian (Eqs.(ld-lf)), which commutes with the S z 
operator. Hence, the eigenvalues for the Izing model can 
be easily calculated without the energy matrix diagonal- 
ization. 

Let us consider the Izing cubic cluster in the case of 
s = 1/2. For each eigenfunction ip m we have found the 
average values (H nn ), (H nnn ) and (H nnnn ) (see Table 1). 

TABLE 1. Matrix elements for the Izing model. 
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Here (H nn ) = {<£m\Hnn\Vm) , {Hnnn) = (<Pm\H nnn \tp m ), 

{Hnnnn) = {fm \H nnnn \fm) , N is the number of the 
wavefunctions (p m with the equal values of {H nn ), {Hnnn) 
and {Hnnnn)- Fortunately, the number of different sets 
of {H nn ), {H n nn), {H nn nn) is much less than the total 
number of the wavefunctions. For example, we have sev- 
enty functions (p rn with S z = and only six different sets 
of corresponding (H nn ), {H nn n), {H nnnn ) (Table 1). 

Using data from Table 1 we can write analitical ex- 
pression for all Izing eigenvalues e m : 

£m = Jl{H nn ) + J2\Hnnn) + Js{H n nnn)- (2) 
Complete set of different e m can be written as follows: 

ef^ = Ji — J% — </3! 

JO) _ T . JO) _ J . JO) T . /o\ 

e 2 — — Ji, e 3 — — j 2 , £4 — —J3, {£) 

e 5 ° = — Ji — J2 + J3; £6°'' = — 3Ji + 3J2 — J3. 
e[ 1] =0.5(Ji- Ji- J 3 ); 
4 1 ' =0.5(-Jx- J 2 + J 3 ); (4) 
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4 1} -0.5(-3Ji+3J 2 - J 3 ). 
4 2) = J i; e< 2) - J 2 ; 4 2) - J 3 . (5) 
ei 3) = 0.5(3Ji + 3J 2 + J 3 ). (6) 

e{ i] =3Ji+3J 2 + J 3 . (7) 

In Eqs.(2-6) the upper index shows S z absolute value, the 
lower index numerates energy levels with the same S z . 

The comparative analysis of Eqs.(3,4) as well as 
Eqs.(6,7) immediately shows that eigenstates with S z — 
±1 or S z = ±3 can not be the ground states It results 
from Eqs. (3,5,7) that the eigenstates with S z = ±2 can 
not be the ground states either. 

Indeed, 4 2) > 4 0) , if Ji > 0. Futher, if Ji < 0, there 

is a competition between 4°^ (i = 1,6) and e^ 4 \ 
Assuming J\ = — 1 and writing the system of inequalities 

' -(0) 



we get 



e?>>-l (i = 1,6); 
4 4) > -1; 



•h £0; 
■h e 0. 



Hence, the energy levels with e\ can not be lowest. 
Analogously, we have from Eqs. (3,5,7) the relations for 

4^ and 4^: 

4 2) > 4 0) , if J 2 > 0, and 4 2) > or 4 2) > e[ 4 \ if 
•h < 0; 

4 2) > 4 0) , if ,h > 0; and 4 2) > or > e[ 4 \ if 
J 3 <0. 

Using Eqs. (3,7) we can find the "phase" boundaries 
between the ground states with S z — and S z = ±4: 

J 3 = - Ji - 2J 2 ; J 3 = -3Ji - 4J 2 ; J 3 = -4J X - 3J 2 ; 



J 3 = -1.5Ji - 1.5J 2 ; J 2 = -Ji; J 3 = -3Ji. (8) 

"Phase" diagram for the Izing cubic cluster with s = 
1/2 is shown in Fig. 2. We use the word "phase" for 
convenience to distinguish between regions with different 
ground states. 

Now we consider the Izing cubic clusters with s > 1/2. 
The analysis for the classical spins (s = oo) is the same 
as for the s = 1/2 since there are only two spin projec- 
tions in both cases. The analitical examination of clusters 
3 5 

with s = 1 ; — , — is simple but cumbersome since there 
' 2' 2 v 

are much more unique sets of parameters (H nn ), (H nnn ), 
{H nnnn ) (66, 129 and 999 for the eigenstates with S z = 
3 5 

and s — 1; - , - , correspondingly). So, we have performed 
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numerical calculations and found that for all values s the 
"phase" diagram is exactly the same as shown in Fig. 2. 

It is interesting to compare the ground state proper- 
ties of the Izing and Heisenbexg cubic clusters. Below the 
case of s = 1/2 is dcscribcd.B For the Heisenberg cluster 
the zero-temperature two-spin correlation functions have 
been calculated using the relation u>ij — (<fto\SiSj\(f>o), 
where </>o is the ground-state eigenfunction. The cubic 
symmetry allows us to take into account only three cor- 
relation functions: loqi, w 2 and loqq. 

First we assume J3 = in Eq.(l). There are four dif- 
ferent cases to be considered: (a) all the interactions are 
ferromagnetic {J\ < 0; J2/J1 = a > 0); (b) all the in- 
teractions are antiferromagnetic (Ji > 0; a > 0); (c) the 
NN interactions are antiferromagnetic, the NNN inter- 
actions are ferromagnetic (J\ > 0; a < 0); (d) the NN 
interactions are ferromagnetic, the NNN interactions are 
antiferromagnetic (Ji < 0; a < 0). 

The case (a) is trivial: the ground-state has 5 = 4, 
all correlation functions are equal to 0.25 irrespective of 
a. In the case (b) for all values of a the ground state 
is non-magnetic (5 = 0), but both signs and magni- 
tudes of the correlation functions are sensitive to the 
a variation (Fig. 3). For example, ujqq changes the sign 
at a = 0.50 and has the highest value 0.25 at a = 1. 
The latter point corresponds to the triplet dimerized 
state with the eigenfunction <fi = (06) (17) (24) (35), where 
(ij) = "/f (Tilj + Tilj)- With a -> +00, UJ02 tends to 
—0.25, ojoi and ujqg approach zero (Fig. 3). 

In the case (c) the ground state is always non-magnetic 
(5 = 0), NN and NNN interactions do not really compete 
(Fig.l), since the opposite direction of NN spins (0 and 
1, 1 and 2, etc) is favorable for the identical direction of 
NNN spins (0 and 2, 1 and 3, etc). As a result, ujq2 is 
always positive 0.25), woi an d W06 are negative and 
depend only slightly on a. 

The case (d) is the most interesting. Fig. 4 shows the 
a-dependence of the correlation functions. If a is small, 
all spins correlate ferromagnetically (loqi = lo^2 = <-^06 — 
0.25). At 0:4 = —0.33 the correlation functions change 
suddenly, so that W02 and W06 become negative and W01 
is lowered by a factor 2.5. It should be noted that \a^\ = 
1/3 with an accuracy of 10 -16 . We found that if | a < 1/3 
the ground state has S = 4. In the case of |a| > 1/3 the 
ground state becomes non-magnetic (5 = 0). 

So far it has been assumed that J3 = 0. If the NNNN 
interactions exist, an analysis becomes more intricate 
since two degrees of freedom appear (a and j3 = J3/J1). 
Fig. 5a shows a "phase" diagram in the case of ferro- 
magnetic NN interactions (Ji = —1) for — 3 < a < 3, 
— 3 < < 3. The line AG separate regions with mag- 
netic (5 = 4) and non-magnetic (5 = 0) ground states. 
The curve AG tends asymptotically to lines J2 = \Ji\ 
and J 3 = 3|Ji| (Fig. 5a). Hence, if either J 2 > |Ji| or 
J3 > 3 1 Ji the ground state is always non-magnetic. Fer- 
romagnetic correlations between spins, which exist in the 
magnetic "phase" in the case of J2 < and J3 < 0, can be 
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destroyed by means of a gradual increasing of Ji and/or 
J3. Let us consider, for example, the case J\ — J2 < 0, 
(3 < 0. If |/3 1 < 1.5, the ground state is magnetic 
(S = 4) and Woi = ^02 = ^06 = 0.25. When |/3| exceeds 
1.5 the ground state becomes non-magnetic (S=0). We 
found that the straight line BC (Fig. 4a), which satisfies 
Ji = J2 = — (2/3) J3 < 0, corresponds to a dimer ground 
state with the eigenfunction <^>o,d = [06] [17] [24] [35] , where 
[km] = 4|(tj4m - Tmlfe)- For the dimer ground state 
the correlation between the NN and NNN spins reduces 
to zero (woi = ^02 = 0), the NNNN spins are coupled an- 
tiferromagnetically (w 6 = —0.75). An another specific 
non-magnetic ground state, corresponding to the line DE 
(Fig. 4a), appears when J\ = J 3 < 0. If a < —0.66 the 
tetramerization takes place. Namely, the non-zero corre- 
lation exists between only the NNN spins: loqi = uqq = 0; 
CJ02 i= 0. The average value of W02 is equal to —0.25 for 
the four-fold degenerate tetramerized ground state (see 
details below). Obviously, the non-magnetic tetramer or 
dimer ground states can be achieved also in a different 
way, simply in limiting cases J2 — > +00 or J3 — > +00. 

The "phase" diagram for J\ = +1 is shown in Fig. 5b. 
The ground state has always 5 = if a > — 1 or /3 > 
—3. The line KM tends asymptotically to lines J2 = 
— Ji and J 3 = — 3 Ji (Fig. 5b) and serves as the boundary 
between magnetic (S=4, loqi = 10^2 = ^06 = 0.25) and 
non-magnetic (S=0) ground states. The dimer ground 
state with the eigenfunction </>o,d occurs in case of a = 1 , 
j3 > 1. The tetramer ground state (woi = ^06 = 0; 
^02 7^ 0) takes place if = 1, a > 1. 

Now we can compare the "phase" diagrams for the 
Heisenberg and Izing models. "Phase" boundaries for 
the cubic s=l/2 Izing cluster are shown in Fig. 5 by dot- 
ted lines. It is important that the Izing "phase" diagram 
can be used as a rough approximation for the Heisenberg 
model (Fig. 5). In particular, the asymptotic straight 
lines for the curves AG and KM (Fig. 5) coincide with 
the Izing "phase" boundaries. The same situation takes 
plase in the case of s = 1.E3 

We have found that the dimer or tetramer ground state 
appears only if a = 1 or /3 = 1, correspondingly. This 
regularity can be explained in the following way. It is 
straightforward to show that H nn = 5*0257 51346 — H nnnn 
and Hnnn = 0.5(5q 25 j + 5i 346 ) - 4sJ> + 1), where 5 02 57 = 

5*0 + 5*2 + 5*5 + 5*7, 5i346 = 5l + 5*3 + 5*4 + 5*6, s = 1/2. 

Since 5 = 5o257 + 5*1346 , we get for a = 1 

H = ,h[0.5S 2 -3 + (/3-l)H nnnn ]. (9) 

We present Eq.(9) for s = 1/2, in the case of s > 1/2 
the numerical parameters should be modified, but the 
structure of the the equation is universal. 

Assuming Ji > 0, it is evident from Eq.(9) that the 
dimerization can be favorable only for (3 > 1, because 
of {(j>o t d\H nnnn \(j)^ d ) = -3 < and 5* = 0. If J\ < 0, 
we have to compare (<f)f\H\(j)f) and (4>Q^d\H\<f>a,d) , where 
4>f = V255, and H is given by Eq.(9). Taking into account 
that (<fif\H nnnn \(j)f) = 1, we have the simple equation 
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(0.5 ■ 20 - 3) + 1 • {0 - 1) = -3 - 3(0 - 1), which results 
to0 = -1.5. 

In the case of = 1, Eq.(l) transforms to 

H = Jl[5 257»S , 1346 + 0.5a(S'o2 5 7 + S 1346 ) - 3a] = 

= MQ.5S 2 + 0.5(a - 1)(S 2 257 + S 2 13m ) - 3a]. (10) 

For all four eigenfunctions <f>t,i (i = 1, 4) of the tetramer- 
ized ground state we have got 

(4>t,i 1^0257 l<At,i> = (0mI S 1346|<? ! O = °- 

Since (</»/|5g 257 |</>/) = (^/l^mel^/) = 6, in the case of 
Ji < the tetramer state has the lowest energy if (0.5 • 
20) + 6 • (a — 1) — 3a < —3a, which leads to a < — |. 
Assuming Ji > 0, it is easily seen from Eq.(10) that the 
tetramer ground state is energetically preferable if a > 1 . 
Otherwise, states with the non-zero averages of S , q 257 and 
S 2 346 are lowest. 

The Hamiltonian in (10) commutes with Sq 257 an d 
S 2 346 . Hence §t,i are the eigenfunctions of these opera- 
tors with the eigenvalues S'0257 = S1346 = 0. To get more 
information about the tetramer ground state we have to 
introduce the dimer spin operator S02 = So + 5 2 with 
eigenfunctions (j) 02 (S 02 , M), which arc defined by the 
equations S$ 2 \<f> 2{S 02 , M)) = S 02 (S 02 + 1)1^02(^02, M)> 
and S 02 , z \(j) 02 {S 02 ,M)) = M\<f) 02 (So 2 , M)) {S 02 = 
0; 1; M = 0; ±1). The similar relations are valid for the 
S57, 5i3 and S^q operators. We have found that 

{4>t,i\ S 01,z\4>t,i) = {4>t.i\ S 57.z\ ( t > t,i) = 
= ((tH,i\Sl3,z\<h,i) = (<Pt,i\S46,z\<f>t,i) =0;i = 1,4 
(hAlfo,!) = (0t,l|^7l^,l> = 

= {4>t.i\Sl 3 \(j)t,i) = {4>t.\\Sl & \cj)ts) = 0; 

((pt,2\ S 02\ < f ) t,2) = (0t,2|S , | 7 |0 t ,2> = 

= (<l>t, 2 \Si 3 \<f>t, 2 ) = ((j>t. 2 \Sl 6 \(f>t, 2 ) = 2; 

(<^,3 1 £021^,3) = (0t,3|S , 57 |^t,3) = 

= (4>t,4\s 13 \<pt,4) = {(l>t,i\s m \(f>t,i) = k ; 

{4>t,i I ^02 1^,4) = (</>i,4|S , 5 7 |^t,4) = 
= (^,3|^3l^,3) = (^,3|S 7 I 6 |^ 3 } = 2-K, 
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where the factor k is fractional and a-dependent. Ex- 
plicit forms for ti 3 and <f>t,4 also depend on the a value. 
The functions 4>t,i and <p t ,2 can be always written as 

0t,l = 002(0, 0)057(0, 0)013(0, 0)^46(0,0); 
0t,2 = ^{002(1, 1)057(1, -1) + 002(1, -1)057(1, 1)- 
002(1,0)057(1,0)} X { -01 3 (1, O)0 46 (1,O) + 
+ 01 3 (1, 1)046(1, -1) + 046(1, 1)013(1, -I)}- 

It is easy to check that (<j)t,i\SoS2\(f>t,i) = —0.75; 
(0t,2|<SoS2|0t,2) = 0.25. Our calculations have showed 

also that always X^i^mI^o^^m) = — 1, and the 
avarage value of u 02 for the tetramer state is equal to 
-0.25. 

In conclusion, we have studied effects of competing ex- 
change interactions on the ground state of the cubic Izing 
1 3 5 

clusters with s=-;l;-;-;oo. The "phase" diagram for 

all Izing clusters is universal. In case of s = 1/2 and s = 1 
we have studied both Izing and Heisenberg clusters. For 
the Heisenberg s = 1/2 cluster it has been found that 
the ground state can be either magnetic with S = 4 or 
non-magnetic with 5 = 0. Intermediate values of the 
total spin S are not realized for the ground state. For 
s = 1/2 the non-magnetic dimer ground state appears in 
two cases: if J\ = J2 < 0; (3 < — | and if J\ = J 2 > 0; 
(3 > 1. The tetramer four-fold degenerate ground state 
takes place if J\ = J 3 < 0; a < — | and J\ = J 3 > 0; 
a > 1 . The two-spin correlation functions are much more 
sensitive to a and (3 changes than the total spin S of the 
ground state. It is important that for s = 1/2 and s = 1 
(and, probably, for higher spin values) the Izing "phase" 
diagram can be used as a rough approximation for the 
Heisenberg model. 
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FIG. 1. A sketch of a cubic cluster with eight spins. 



FIG. 2. The Izing cubic cluster. Dotted lines 

show "phase" boundaries between the ground states with 
S z ,max = ±8s (±4 for s = 1/2, ±8 for s = 1, etc.) and 
S z =0. 

FIG. 3. The Heisenberg cluster with s = 1/2. The 
a-dependence of the two-spin correlation functions: Ji > 0; 
J 3 = 0. 

FIG. 4. The Heisenberg cluster with s = 1/2. The 
a-dependence of the two-spin correlation functions: Ji < 0; 
J 3 = 0. 

FIG. 5. The Heisenberg cluster with s = 1/2. "Phase" 
diagrams for (a) Ji < 0; — 3 < a < 3; —3.5 < j3 < 3, and 
(b) Ji > 0; -3.5 < a < 2; -5.5 < (3 < 2. Dotted lines 
show boundaries between the ground states with S z = ±4 
and S z = for the Izing model. 
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